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Abstract 

The search for charged Higgs bosons, which are predicted in supersymmetric 
theories, is difficult at hadron colliders if the mass is large. In this paper we present 
the theoretical set-up for the production of charged Higgs boson pairs at the LHC 
in gluon-gluon collisions: pp— > H + H~ . When established experimentally, the 
trilinear couplings between charged and neutral CP-even Higgs bosons, H + H~h° 
and H + H~ H , can be measured. 



CERN-TH/97-137 
DESY 96-254 
[hep-ph/9707430 
July 1997 



1 Introduction 



Charged Higgs bosons are part of the extended Higgs sector predicted in supersymmetric 
theories. They belong to the two Higgs isodoublets that must be introduced to provide 
masses to down- and up-type fermions through the superpotential and to keep the theory 
free of anomalies JT|. The mass of these particles is expected to be in the range of the 
electroweak symmetry-breaking scale v = 246 GeV, albeit with a rather wide spread. 
Embedding low-energy supersymmetry into supergravity models with universal SUSY 
breaking, the masses of the heavy neutral and charged Higgs bosons are indeed predicted 
at the TeV scale for large parts of the parameter space [see e.g. Ref. 0]. 

At hadron colliders ||, several mechanisms can be exploited to search for charged 
Higgs bosons. A copious source of charged Higgs bosons with fairly light mass are decays 
of top quarks, t — > b + H + . Since top quarks are produced with very large rates at the 
LHC, charged Higgs bosons can be searched for in this channel for masses up to 

the top mass. No experimental technique has been established so far for charged Higgs 
bosons with masses above the top mass.[] Channels such as pp — > ibH + , in which the 
Higgs bosons are emitted from heavy quark lines and the Drell-Yan production of 
charged pairs ||, pp — > qq — > H + H~ , are potential candidates for the search. 

Extending the previous analysis of neutral Higgs bosons in supersymmetric theories, 
Ref. ||, we present in this paper the theoretical set-up for the production of charged Higgs 
boson pairs in gluon-gluon collisions: 

PP~> 99 -> H + H~ . (1) 

The analysis of this channel is motivated by the large number of gluons in high-energy 
proton beams. We work out the predictions for the cross sections^ at the LHC. For 
charged Higgs masses above the top mass, yet less than 300 GeV, the cross sections will 
turn out to be in the range between 1 fb and 10 fb if the mixing parameter tan/? in 
the 2-doublet Higgs sector is either small, tan/5 > 1, or large, tan/? ~ 50. The search 

1 At e + e~ colliders, charged Higgs bosons can be searched for in the process e + e~ — * H + H~ for masses 
up to the beam energy, i.e. up to m H ± < ITeV for prospective 2 TeV e + e~ linear colliders. 

2 The triangle contributions to the production cross section have already been determined in Ref. 0, 
but the more complicated box contributions were only estimated in this paper. 

3 After finalizing this calculation we received a copy of the paper The present analysis differs from 
this paper in several aspects: (i) We present the partonic cross sections in compact analytical form, while 
they are given in Ref. JsJ in terms of twelve form factors from which the individual UV singularities are 
not separated yet. (ii) We include a careful analysis of the theoretically interesting large quark-mass 
limit, which serves as an important cross-check of the calculation, (iii) Our final results differ from those 
in Ref. ||: Part of the discrepancies could be traced back to a wrong flux factor and a wrong factor 
of 2 in the gluon luminosity of || (moreover, opposite to the calculation in Ref. ||, Z-exchange cannot 
contribute to the gg fusion to charged Higgs boson pairs as a result of symmetry arguments) ; other points 
could not be isolated in detail since the cross sections in [|| are evaluated numerically in a global way 
that does not allow checks of separate terms. The results presented in our paper have been cross-checked 
in two independent calculations. 
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for such events will therefore be difficult, and the characteristic decay properties must 
be exploited exhaustively in the experimental analyses^. However, if the signal can be 
established experimentally, the trilinear couplings between charged and neutral CP-even 
Higgs bosons, h° H + H~ and H° H + H~ , can be studied, thus complementing related 
analyses in the purely neutral sector || |J. Measurements of these couplings are the first 
step in the reconstruction of the self- interaction terms in the Higgs potential. By contrast, 
the emission of charged Higgs bosons from top quarks and the Drell-Yan production of 
charged Higgs pairs do not include the trilinear couplings. 

The paper is organized as follows. After a brief summary of the charged Higgs sec- 
tor within the Minimal Supersymmetric Standard Model (MSSM) in the next section, 
the theoretical set-up for the production of charged Higgs boson pairs in gluon-gluon 
collisions at the LHC will be described in the subsequent section which includes a brief 
phenomenological evaluation. 



2 Charged Higgs Bosons: SUSY Expose 



The MSSM will be adopted as the paradigm for a theory of charged Higgs bosons. After 
the characteristic phenomena are worked out for this example, the results can easily 
be adapted to more general scenarios. At least two Higgs isodoublets are required in 
supersymmetric theories to generate masses for down- and up-type fermions through 
the trilinear superpotential, and to cancel anomalies associated with higgsino triangle 
loops. After absorbing three of the eight field components to generate the longitudinal 
components of the electroweak gauge bosons, three neutral CP-even and CP-odd particles 
and a pair of charged particles build up the Higgs spectrum of the MSSM 

In the minimal version, the properties of the Higgs particles depend essentially on 
three parameters: two masses, generally identified with the Z-boson mass mz and the 
pseudoscalar Higgs-boson mass rriA, and the ratio of the two vacuum expectation values of 
the neutral Higgs fields, tan/5 = v 2 /vi, assumed to vary between unity and m t /m b ~ 50. 
Including the leading radiative corrections, induced in the potential by the large top mass, 
the other masses can be expressed in terms of mz and m^. 



m h°,H° ~ o 



m AZ ~F 



ni 



AZ 



Am\m 2 z cos 2 {2(3) — 4e (m\ sin 2 /? + m 2 z cos 2 /3) 



(2) 
(3) 

with the abbreviation m\ z = m\ + m 2 z + e ; the leading corrections are characterized by 
the radiative parameter e [|10|j , 



m 2 H ± 



2 i 2 
m A + m W 



3G f 



V2tt 2 sin 2 /3 



log 



m 



771 



(4) 



4 Such an experimental analysis is beyond the scope of the present paper, which is restricted to the 
theoretical basis of the gluon-fusion process. 
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The parameter m~ = m^m^ denotes the average mass square of the stop particles. The 
mixing in the CP-even Higgs sector is given by the angle a with 



m A + m 



tan (2a) = — £ f tan (2(3) -- < a < . (5) 



m\ — m| + e/cos (2/3) 



7T 



2 



The mass pattern is a consequence of the fact that the size of the quartic self-couplings 
is determined by the square of the gauge couplings in the MSSM. Since these are small, 
strong upper bounds on the mass of the lightest neutral CP-even Higgs boson h° can be 
derived. The lightest neutral Higgs boson h° may therefore be the first Higgs particle to 
be discovered. However, since it is very difficult to determine tan/3, the measurement of 
the h° mass will in general not enable us to predict the mass of the charged Higgs bosons. 
Moreover, as a result of the decoupling theorem, m^o becomes increasingly insensitive to 
the values of the heavy Higgs masses beyond ~ 200 GeV, as is evident from Fig. 1. 




60 100 200 400 

m A [GeV] 



Figure 1: The masses of the Higgs bosons in the MSSM for two representative values of 
tanf3 =1.5 and 30. 



The couplings of the SUSY Higgs particles to fermions are modified by the mixing 
angles with respect to the Standard Model (Table 1). They are renormalized only indi- 
rectly in leading order through the renormalization of the mixing angle a. By contrast, 
the trilinear Higgs couplings are modified indirectly through the renormalization of the 
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SM 


9t 
9b 


(V2G F y 2 m t 
(V2G F y 2 m b 


MSSM 


tth° 
bbh° 
itH° 

bbH° 
ibH + 
btH- 


(cosa/sin/3) g t 
— (sina/cos/3) g b 

(sina/sin/3) g t 

(cosa/cos/3) g b 
-± [g t cot/3 (1 - 75) + 9b tan/3 (1 + 75)] 
-j- [g t cot/3 (1 + 75) + g b tan/3 (1 - 75)] 



Table 1: The Higgs- quark couplings to fermions in the MSSM relative to the corresponding 
SM couplings. 



mixing angle a, but also explicitly through vertex corrections [II]. Normalized in units 
of Ao = (\f2Gp) 2 Tn%, they are given by 

X(h°H + H-) = 2 cos 2 e w sin {(3 - a) + cos (2/3) sin (J3 + a) + * ^^f^ (6) 



\{H°H + H-) = 2 cos 2 9 w cos (/3 — a) — cos (2/3) cos (f3 + a) + 



Mf sin/3 

e sina cos 2 /3 
M| sm^ 



(7) 



Their size is illustrated for tan/3 = 1.5 and 30 in Fig. 2. The trilinear couplings of the 
light scalar Higgs boson to charged Higgs particles turn out to be positive in the relevant 
charged-Higgs mass range, while those of the heavy scalar Higgs boson are negative for 
small charged Higgs masses and change to positive but small values for charged Higgs 
masses between about 150 and 250 GeV. The absolute size of these couplings extends 
from (9(10 -1 ) to 0(1). For large values of tan/3 the trilinear couplings change rapidly 
near m#± ~ 150 GeV. 

When the form factors are evaluated numerically, the masses and couplings are con- 



sistently used in two-loop order [12 



3 Charged Higgs Pairs in gg collisions 

Two mechanisms contribute to the production of charged Higgs-boson pairs in gg fusion, 
exemplified by the generic diagrams in Fig. 3: (i) Virtual neutral CP-even Higgs bosons 
h°, H°, which subsequently decay into H + H~ final states, are coupled to gluons by t, b 
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Figure 2: The trilinear couplings of the charged Higgs bosons in the MSSM for two repre- 
sentative values of tanf3 = 1.5 and 30. 



quark triangles; (ii) The coupling between charged Higgs bosons and gluons is also me- 
diated by heavy-quark box diagrams. As a result of CP invariance, the neutral CP-odd 
Higgs boson A does not couple to pairs of charged Higgs bosons. The Z boson cannot 
mediate the coupling either: The vector component of the wave function does not couple 
to the initial gg state as a result of the Landau- Yang theorem, the CP-odd scalar com- 
ponent does not couple to the CP-even H + H~ state. Forbidden by the Landau- Yang 
theorem, virtual photons cannot contribute either. 

In the triangle diagrams of Fig. 3 the gluons are coupled to the spin S z — along 
the collision axis. The transition matrix element associated with this mechanism can 
therefore be expressed by the product of one form factor F®, depending on the scaling 
variable tq = Attiq/s, with the generalized coupling defined as 

C*A = ^H+H-Hi~ 9 ~- fi 9q* ■ (8) 

^ s-mir +im H .T H . ^ 

The couplings g^ denote the Higgs-quark couplings in units of the SM Yukawa couplings, 
collected in Table 1; the Higgs self-couplings \n+H-Hi have been defined in Eqs. (6), (7). 
VI is the gg center-of-mass energy. For the sake of convenience, the well-known triangle 
form factor F^(tq) is recapitulated in Appendix 1. In the limit of large [t] and small [b] 
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Figure 3: The s-channel triangle diagrams and box diagrams contributing to gg — > H + H . 

quark- loop masses, the form factor simplifies considerably (r t = 4m 2 /s, r b = 4m 2 /s): 

il(rt»l) - | (9) 

(10) 

These expressions provide useful approximations in practice for a large range of Higgs 
masses. 

The box diagrams in Fig. 3 contribute to the two spin states S z = and 2. The tensor 
basis for the states is given explicitly in Appendix 1. For S z — 2 two independent elements 
can be defined, carrying positive and negative parity under space reflection. For 5*2 = 
only the basis element carrying positive parity is realized since the [S z = 0; P = — ] 
element is odd under CP transformation and thus forbidden for CP conserving gg — > 
H + H~ transitions. Therefore we are left with three independent form factors F, G and 
if, corresponding to the spin-C/P states ++ , 2 ++ and 2 . It turns out to be convenient 
to split each form factor into three components, which are proportional to tan 2 /? and 
cot 2 /?, and independent of tan/?: 

F n = tan 2 /?F^ + cot 2 /?F^ + F^ (11) 
G u = tan 2 /?G^ + cot 2 f3G 2 D + G 3 D (12) 
H u = tan 2 /?#i + cot 2 13 + . (13) 

The form factors depend on the square of the invariant energy s and the momentum 
transfer squared t, u. These Mandelstam variables are given by the gluon beam energy E 



F h A {n « l) 



n 



log — + in 
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and the H ± production angle 9 with respect to the gg axis in the partonic cm. system: 



s = AE 2 

t = m 2 H+ - 2E 2 (l - $ cos#) ( 14) 

u = m 2 H+ - IE 2 (l+$ cosfl) 

with the velocity (3 = (1 — 4m 2 H± /s) 1 / 2 . Compact analytical expressions in terms of these 
Mandelstam variables are given in Appendix 1. 

In the double limit of large t mass and small b mass, the form factors can be reduced 
to very simple expressions: 



*s - 


-> 






F 2 u - 


n 

4 


log — + 111 

L 4 


2 

+ 


F D - 


n 
2 


log — + 27T 


2 


G n - 


-> 






H n - 


-> 0. 







As expected, the S z — 2 form factors Gn and ifn vanish in the limit of large top 
masses. The non-zero components of the S z — form factor Fo can be interpreted in 
a transparent way. Pinching the top propagator in the bbbt box effectively leads to a 
bottom-quark triangle, which is accounted for by the r-dependent expressions in Fo- The 
remaining constant term 2/3 corresponds to the sum of the asymptotic limits of the bttt 
box [= —1/3] and the contribution generated by the bbtt box [= +1]. 

The form factors are compared with their asymptotic values in Fig. 4 for a particular 
value of the scattering angle. The partonic center-of-mass energy has been chosen near 
threshold, = 310 GeV for m H ± = 150 GeV, where the dominant contributions to the 
hadron cross section are generated. The imaginary part of Fa is very small and is therefore 
not shown in the figure. Apparently, the box form factors approach the asymptotic values 
only for large quark masses. For m t — 175 GeV the asymptotic box expressions are not 
yet useful in practice. Notably the real part of Fo, which for most of the parameter 
space gives the leading contribution to the cross section, approaches asymptotia only 
very slowly. Nevertheless, the comparison in the asymptotic parameter range provides an 
important cross check of the general calculation for realistic b, t quark masses. 
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The differential cross section at the parton level can finally be written in the form 



f b» -^1 = ^1 [I Ecg*2+*b|' + IGbl" + !*□!'] • us) 

The hadronic cross section for the process pp —>■ gg —>■ H + H~ can finally be derived by 
integrating ( |TED over the production angle < 9 < n and the gg luminosity: 

a [ pp ^ gg ^ h+h-] = / dr— a(s = rs) . (17) 

JAm 2 H± /s " r 

It is likely that QCD corrections to gluon-fusion processes, known for the triangle di- 
agrams, unknown still for the box diagrams, enhance the event rate (cf. Ref.jnj). In 
parallel to other gluonic processes, the corrections are expected to remain under control 
if the factorization scale in the gg luminosity is chosen to be of the order of the invariant 
(H + H~) energy. 

We have evaluated the pp cross section numerically^ for the LHC energy y/s = 14 TeV. 
The (pole) mass of the top quark has been set to m t = 175 GeV and the bottom quark 
mass to nib = 5 GeV. The results are shown for three values tan/3 = 1.5, 6 and 30 in Fig. 5 
as a function of the charged Higgs boson mass. [The masses of the neutral Higgs bosons 
h° and H° which enter through the charges C^, are fixed by m H ± and tan/3 in the MSSM.] 
For small values of tan/3 and modest values of m H ±, the triangle diagram provides the 
dominant contribution. With rising Higgs mass and for large tan/3, the box diagrams 
become increasingly important; they dominate for Higgs masses above the top mass. The 
choice tan/3 = 6 in Fig. 5 corresponds to the minimal value of the cross section. This is 
demonstrated in Fig. 6 for a fixed Higgs mass m#± = 200 GeV and tan/3 varied between 1 
and 50. The minimal cross section is realized at tan/3 ~ yjnit/mb ~ 6, roughly equivalent 
to the minimum of the average QQH ± coupling (gQQH±) = { m t tan 2 /3 + m 2 cot 2 /3)^ 2 . 

The cross section for gluon fusion of charged Higgs boson pairs is smaller than the cross 
section for the Drell-Yan mechanism. In the interesting mass range m^± > 200 GeV, 
the ratio is about 0.4 and 0.2 for tan/3 = 30 and 1.5 respectively. It is noteworthy that 
the angular distributions of the charged Higgs bosons are different for the two mech- 
anisms. While the Drell-Yan mechanism generates a sin 2 # distribution in the (H + H~) 
center-of-mass frame, the gluon-fusion processes are not suppressed by angular-momentum 
conservation in the forward direction; for example, the amplitudes built up by the Higgs 
exchange diagrams, are isotropic in 9. Even though the gluon-fusion mechanism is not 
dominant for the production of charged Higgs-boson pairs, it is nevertheless of great in- 
terest since the gg fusion cross section and the angular distribution are affected by the 
size of the trilinear self-couplings between the CP-even neutral and charged Higgs bosons. 



5 FORTRAN code in Rcf. F3 
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100 150 200 250 300 350 400 

m H ±[GeV] 

Figure 5: Total cross sections (full lines) and the contributions without the triangle dia- 
grams (dashed lines) for gluon fusion of charged Higgs-boson pairs at the LHC [CTEQ4M 
parton densities I^Bjj with a s (m 2 z ) = 0.116]. 



11 




1 2 3 4 5 6 7 8 10 20 30 40 

tan(3 



Figure 6: Dependence of the gluon-fusion cross section of charged Higgs-boson pairs on 
the parameter tanf3. 

4 Conclusions 

As shown in the representative examples of Fig. 5, the cross sections for charged Higgs 
masses below 300 GeV vary between 1 and 10 fb in a large part of the MSSM parameter 
space. This corresponds to rates of 100 and 1,000 events for an integrated luminosity 
of J C = 100 fb -1 , which is expected in a year of high-luminosity runs at the LHC. For 
small Higgs masses the triangle contribution appears large enough to generate a significant 
transition amplitude; this will enable us to measure the trilinear self-couplings Xn+H-h° 
and \h+h-h° of the CP-even neutral and charged Higgs bosons. Charged Higgs parti- 
cles decay predominantly to rv T final states at large tan/3, leading to asymmetric t + t~ 
configurations in the decays of the two Higgs particles. For small tan/3, the potentially 
dominant chargino/neutralino decays are more difficult to control experimentally. These 
experimental problems must be approached, however, in dedicated detector and back- 
ground simulations, which are beyond the scope of this first-step theoretical analysis. 
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APPENDICES 



1. Form factors in gg — ► H + H 
Parameter definitions : 

S = (Pa+Pb) 2 , i=(Pa-Pc) 2 , U=(p b -p c ) 2 

ti — t — m 2 H ±, ui=u — m 2 H ±, p 2 — t\U\/s — m H ± 
tq = 4m 2 Q /s, $=(1- H±/s) 1/2 

Scalar integrals 



a/S'y 



d 4 q 



D 



ct/37<5 



in 2 [qi - m 2] [( q + p .) _ m 2] [( q + p . + pj j _ m 2] 
d A q 1 



[q 2 - ml] [(q + p^ - m|] [(g + Pi+Pj) - m^] [(g + p, + p d + p h ) -m 2 ] 



with a,P,j,6 — t,b and z, j, k = 1, . . . , 4 

Pi = Pa , P2 = Pb , P3 = -Pc , Pa = ~Pd in a + 6 -> c + <i 



Tensor basis: 



V i 1 

P = + S z = : Af = ^ P " Pb 



(PaPb) 

C _ O 4^ _ , PiPaPb 2 (PbPc) PgP^c 2 ( PaPc ) jffig 2/£# 

— Z . S±2 — g ~\ 2~7 \ 27 \ 27 \ 



P = - & = 



PtiPaPb) P 2 t {PaPb) PtiPaPb) Pt 

Af = t -e w « : CP odd not coupling to H + H- 

\PaPb) 

S - 2 ■ - Pc eUPaPbPc + p v c e» PaPbPc + (p b p c ) ^ vpaPc + {p a p c ) e^ PbPc 

4 {PaPb)p 2 t 

normalization and orthogonality: AiAj = 
Matrix element : 

M (g a9b -> tf+/r ) = .Ma + M.n 

M * = ^rH ClFlA x ^ ah 

2 ^ 27V Q=t,b 
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Triangle form factor: 



F2 = T Q [l + (l-T Q )f(r Q )} 



f(r Q ) 



-i 2 



arcsm 



2 1 



TQ<1 



Box form factors: 



F, 



2ml 



□ 



(1 + 2m? Cg) + (mf + mg - < ± ) { - |c 



*1 /-«tt6 u l /-rttft 



13 - ^23 + m t \ U 12Z 



b <-> t 



= Fq (m t «- m b ) 



= 2m 2 b m 2 t 



2 2 

+ 2m b m t 



4tti 



< ri' - -Mf + - 1) (n™ + ds*) + {f K + ™* + - 1} 



6 t 



m^ 
p\ L 



+ m 2 -m 2 H± ^(p 2 



M 



2 

+{M 6 + L bt (t\,£ 1 )}£f 2 f + {M 6 + ^ ^ " r 



+ 



+ (m 2 + m 2 )p 2 (\-^f)}D\^ 



+ 



mi, 
p\ L 



b <-> t 



where: M bt = ^ (m? — m 2 H ±) 2 — m\ +2(m 2 — m 2 H ±) 



M 6 = m ° ~ ~ + < + m 2 m 2 (2*^1 - l) 



s L S s 



L bt (ii,ui) = m\ 



1 Ui 

2 + T 



0-f) 



S J s 



, 4 

+ m t 



1 2u + m^± 

2 I 



+ m h 



P 2 Pt m H± 



+ m 2 



u 



r (2« + t) - 



Pt 



u 
s 



+ u m H± 
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G a = G D (m t <-> m b ) 



3 _ 4m t 2 m 2 . 



% 2 



- (2m 2 - 2m 2 . + 2m #± ~ ») C 



ttt 

12 



2tit 



2ni u 



tbt 
34 



+ L bt (u)D\f 3 + L w (t)£>£* + { {ml - m\f + (m t 2 + m 2 ) p t 2 }z) 



tt66 
132 



+ 



4m 2 m 2 
sp 2 



6 t 



where: L bt (u) = (m 2 — m 2 ) 2 + -u (u — 2m 2 ) — 2m 2 -u 2 /s and 



□ 



m 2 r 
P? L 



(t! - fix) M 2 Cft + § (t + m| ± ) Cf 3 & - ^ (« + < ± ) C§> + 



" 1_tl V 2 ctf 



+ *«(fi, U!)Dg« - + (^i) { (m 2 - m 2 ) 2 + (m 2 + m 2 ) p 2 }l) 



tt&6 
132 



+ 



m^ 
P 2 t L 



where: R bt {h,ui) = - (u + m H ±) — + - 



UUi 1 (u\—t\ 



m t — m. 



2\ 2 



rn 2 



+ m 2 



" /i - \ ft 



M, 



2 m H± + rnf — ml 1 



6t 



2 

tj-2 _ m i rrl 

H n — 2 

m 2 



2. Pinching top propagators 



Because of delicate cancellations, the limit m t — > oo is not easy to calculate for the box 
diagrams involving one t line and three b lines, etc. We have performed this limit for the 
amplitudes first in the Feynman parametrization. The problem however can also be solved 
in a systematic way by adopting the scheme for the expansion in the heavy-quark mass, 
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which has been discussed in Ref. [[HJ. This method is technically easier than the Feynman 
parameter method. The limit is performed by expanding the scalar integrals inside the 
form factors; the m t contributions arising from the Yukawa couplings are factored out. 



The procedure |17| will be illustrated for a scalar three-point function including one 
top propagator: 



(-a 



bbt 

felfc2 



1 



in 2 [q 2 



m; 



;] [(? + h 



m. 



;] [(q + k 1 + k 2 ) 2 - ml] 



which must be expanded up to 0((kikj) 2 jva\). The integral is apparently ultraviolet and 
infrared finite; nevertheless, we will work in n dimensions in order to properly define the 
divergences that occur at intermediate steps of the calculation. Since the Taylor expansion 
with respect to the heavy quark mass 



1 



(q + h + k 2 ) 



(q + h + k 2 ) 



+. 



mi 



mr 



ml 



Nr 



leads to ultraviolet divergences of the integral at q — > oo, the expansion cannot be applied 
in this naive form. However, after regularizing the divergences in n dimensions, Nyy can 
be used in the identity 



1 



(q + ki + k 2 ) - m\ 



N uv + 



1 



(q + ki + k 2 



m 2 t 



-N 



uv 



RlR 



to define a kernel 



R 



IR 



(q + h + k 2 ) 



[(q + h + k 2 ) 2 



which, together with the two bottom-quark propagators in C, can be expanded in the 



variables (qh) / (q 2 



m. 



(kikj)/m 2 and m\jm 2 t . After performing the integral of the 



expansion in n dimensions, the leading terms generate singularities which exactly cancel 
the singularities of the integral over N uv . The finite pieces and the subsequent terms of 
the expansion lead to a well-behaved series in inverse powers of the heavy quark mass. 
The expansion for this particular example can finally be written: 



(~ibbt 



m% 
1 



m 



+ o 



t L 



[kikj / 



mt 



mi 



l + log^ij-^) 
mr. I 1 + 2 log ( ^| ) -g(n b ) 



) + | 
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where = 4m 2 I ^1 anc l 



log T 



Z7T 



2 — 2-y/Kb — 1 arcsin —L 



: « 6 < 1 
: K b > 1 



The same procedure can be applied to the large quark mass limit of scalar three- 
point functions including two top propagators. The product of the two top propagators 
is first naively Taylor expanded. The top part is then split into UV- and IR-divergent 
contributions as shown before. After the Taylor-expansion of the IR-divergent part, the 
calculation of the integrals leads to: 



rittb 



mi 



ml 



m. 



i+ (4)1 



3 2 2 



+ o 



(kikj) 



mi 



For the scalar three-point function involving three top propagators, the naive expan- 
sion is UV-finite, so that it need not be split. One obtains the following expansion up to 
the required order 0((kikj) 3 /m^): 



1 1 

2mf 



m. 



kiki + k\k 2 + k 2 k 2 



12 



mi 



(k 2 + k^hh k\k 2 



30 



+ 



36 



+ 



(hh) 2 

45 



+ 



60 



+ o 



rn) 



The expansion of the four-point functions can be carried out in an analogous way. 
The result is still too long to be presented in plain text; yet the analytic expressions are 



available from 14 



In contrast to the expansion for large t mass, the subsequent expansion for small b 
mass m 2 /s -C 1 is straightforward. 
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